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Abstract 



Oh! 

. We present a new method for determining the renormalization of Green functions to all orders 

- ^ , in perturbation theory, which we call the displacement operator formalism, or the D-formalism, in 

X 

^ I short. This formalism exploits the fact that the renormalized Green functions may be calculated 

by displacing by an infinite amount the renormalized fields and parameters of the theory with 
respect to the unrenormalized ones. With the help of this formalism, we are able to obtain the 
precise form of the deformations induced to the Nielsen identities after renormalization, and thus 
derive the exact dependence of the renormalized Green functions on the renormalized gauge-fixing 
parameter to all orders. As a particular non-trivial example, we calculate the gauge-dependence 
of tan (j at two loops in the framework of an Abelian Higgs model, using a gauge-fixing scheme 
that preserves the Higgs-boson low-energy theorem for off-shell Green functions. Various possible 
applications and future directions are briefly discussed. 

PACS numbers: ll.15.Bt. ll.lO.Gh 
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I. INTRODUCTION 



Renormalization is of central importance in quantum field theory [1-12] , and its consistent 
implementation has been crucial for the advent and success of gauge theories in general, and 
of the Standard Model of strong and electroweak interactions in particular. Despite half a 
century of practice, however, renormalization remains a delicate procedure, mainly because it 
interferes non-trivially with the fundamental symmetries encoded in the Lagrangian defining 
the theory. The subtleties involved manifest themselves at almost every step, ranging from 
the necessity to employ regularization methods respecting all relevant symmetries, the need 
for renormalization schemes that do not spoil the important constraints imposed by the 
symmetries on the Green functions of the theory, all the way to the practical, book-keeping 
challenges appearing when the renormalization is carried out in higher order calculations. 

In this paper we develop a new formalism, which we call the Displacement Operator 

Formalism, or the D-formalism, in short. The D-formalism enables one to systematically 

organize and explicitly compute the counterterms (CTs) involved in the renormalization 

procedure, to all orders in perturbation theory. The central observation, leading to this new 

formulation, is that the effect of renormalizing any given Green function may be expressed 

in terms of ultraviolet (UV) infinite displacements caused by the renormalization on both 

the fields and the parameters of the theory. Specifically, these UV infinite displacements, or 

shifts, quantify the difference between fields and parameters before and after renormalization. 

For example, in the case of a theory with a scalar field (f), and two parameters, the coupling 

constant A and the squared mass m^, the corresponding UV infinite shifts S(f), 5X, and 

1 

are given by 5(f) — {Z^ — 1)0r, SX — {Z\ — 1)Ar, and Srn^ — {Z^i — 1) m|j, where 

1 

the renormalization constants are defined, as usual, through = ^^^r, A = ^aAr, and 
= Zm2 rn^. Evidently, in this formulation, dynamical fields and parameters are treated 
on a completely equal footing. 

In order to systematically expose the way in which these shifts implement the renor- 
malization at the level of Green functions, one introduces the displacement operator D, a 
differential operator given hy D — S(f>-^ + ^"^^^ + '^"^^a^- turns out (see section II for 
details) that the net renormalization effect is captured to all orders by the exponentiation 
of the D operator. Thus, the renormalized (n-point) Green functions are eventually ob- 
tained from the bare ones, T^n, through the master equation (p^T^n = (e^(0^r<^n)), where 
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the brackets (...) means that the corresponding shifts, imphcit in the D operator, are to be 
replaced by their expressions in terms of the renormahzation constants given above, after 
the end of the differentiation procedure. 

It is important to appreciate at this point the inherent non-perturbative nature of the 
above formulation, manifesting itself through the exponentiation of the D operator. Per- 
turbative results (at arbitrary order) may be recovered as a special case through an ap- 
propriate order-by-order expansion of the above master formula, whose validity however is 
not restricted to the confines of perturbation theory. This is to be contrasted with other 
well-known renormahzation methods (as, for example, the "algebraic approach" [13-15]), 
which are formulated at the perturbative level. 

One of the main advantages of this formulation is the ability it offers in determining 
unambiguously the CTs to any given order in perturbation theory. Specifically, the con- 
tributions of the CTs are automatically obtained through the straightforward application 
of the D operator on the unrenormalized Green functions, without having to resort to any 
additional arguments whatsoever. This last point is best appreciated in the context of com- 
plicated non-Abelian gauge theories, or special gauge-fixing schemes, where keeping track of 
the CTs can be not only logistically demanding, but at times also conceptually subtle. For 
example, possible ambiguities related to the gauge-fixing parameter (GFP) renormahzation, 
endemic in sophisticated quantization schemes such as the Background Field Method (BFM) 
[16-20] are automatically resolved (see section III). As we will see in detail in the main body 
of the paper, the fact that the D-formalism correctly incorporates the action of the CTs is 
clearly reflected in the absence of overlapping divergences from the resulting expressions. 

An important property of the D-formalism is that it reproduces exactly the usual di- 
agrammatic representation of the CTs, if one acts with the D operator on the Feynman 
graphs determining the given Green function before the integration over the virtual loop 
momenta is carried out. If one instead acts with D after the momentum integration has 
been performed one loses this direct diagrammatic interpretation, but recovers the same final 
answer for the renormalized Green function. This last point becomes particularly relevant 
in the case of gauge theories, where one of the parameters that undergoes renormahzation 
is the GFP, to be denoted by ^. This, in turn, will introduce in the D operator a term of 
the form 3^-^, which must act on the corresponding Green function. Clearly, for this to 
become possible the dependence of the Green function on ^ must be kept arbitrary, that 
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is, one may not choose a convenient value for ^, like for example ^ = 1. Evidently, if one 
were to first carry out the momentum integration and then differentiate, one would be faced 
with the book-keeping complications of computing with an arbitrary ^. If, instead, one opts 
for the action of D on the Green function before the loop integration, the straightforward 
use of the chain rule for the strings of tree-level propagators (and possibly vertices), and a 
subsequent choice of, say, ^ = 1, is completely equivalent to the conventional approach of 
computing at a fixed gauge. Obviously, one may choose either one of the two procedures, 
i.e., integrate before or after the action of depending on the specific problem at hand, a 
fact which exemplifies the versatility of the new formalism. 

Perhaps one of the most powerful features of the D-formalism consists in the fact that it 
provides a handle on the structure and organization of the CTs to all orders in perturbation 
theory. In particular, it furnishes the exact form of the deformations induced due to the 
renormalization procedure to any type of relations or constraints which are valid at the 
level of unrenormalized Green functions. Such deformations will arise in general if the 
symmetries (global or local) enforcing the relations in question involve static parameters, 
which cannot admit kinetic terms compatible with those symmetries. A well-known example 
of such a type of relation is the so-called Nielsen identities (NI) [21], which, by virtue of 
the extended Becchi-Rouet-Stora (eBRS) symmetry [22, 23], control the GFP-dependence 
of bare Green function in a completely algebraic way. The NI are very useful in unveilling in 
their full generality the patterns of gauge cancellations taking place inside gauge independent 
quantities, such as ^'-matrix elements [12, 24]. Gontrary to the case of the usual BRS 
symmetry [25], however, where only dynamical fields are involved, the eBRS symmetry 
involves also static parameters (^, r/) which cannot be promoted to dynamical variables 
without violating this latter symmetry. Consequently, the NI are affected non-trivially by 
the ensuing renormalization, which deforms them in a complicated way [26] . As we will see 
in detail in section IV, the straightforward application of the D-formalism yields, for the 
first time, the deformation of the NI in a closed, and, in principle, calculable form. This is a 
clear improvement over the existing attempts in the literature, where the deformations are 
formally accounted for, and their generic structure and properties inferred, but, to the best 
of our knowledge, no well-defined operational procedure for their systematic computation 
has been spelled out thus far. 

The paper is organized as follows: In section II we introduce the D-formalism in the 
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context of a scalar (p^ theory, and explain its most characteristic features. The generalization 
of this formalism to more complicated field theories is presented, and a general master 
equation valid for any Green function is derived. In section III we apply the £)-formalism to 
the case of known examples, in order to familiarize the reader with its use. In particular, we 
study the two-loop rcnormalization for the cases of (j)'^, QED, and QCD formulated in the 
BFM, and demonstrate how all necessary CTs are unambiguously generated through the 
straightforward apphcation of the D operator. Section IV contains the first highly non-trivial 
application of the D-formalism, namely the derivation of the deformation induced to the NI 
by the rcnormalization procedure. First, we review the derivation of the NI, and introduce a 
consistent graphical representation. Then, by means of the D-formalism, we obtain for the 
first time the NI deformation in a precise, closed form, presented in Equation (4.25), which 
is a central result of this paper. 

In section V we revisit the Abelian Higgs model [27], and study its quantization in the 
special class of gauges which preserve the Higgs-boson Low Energy Theorem (HLET) [28- 
32] beyond the tree-level. The reason for turning to this model is three-fold. First, it has a 
rich structure, most notably two different GFPs, thus providing an ideal testing ground for 
the newly introduced formalism. Second, in theories with spontaneous symmetry breaking 
(SSB), the HLET is of central importance, and together with the Slavnov- Taylor identities 
[33] and NI, severely constrains the form of the Green functions of the theory. By quantizing 
the model using the HLET- preserving gauges [34] we allow for the consistent implementation 
of a very special rcnormalization condition for the Higgs vacuum expectation value (VEV) 
V. Specifically, if the gauge fixing does not preserve the HLET, as is, for example, the case 
of the ordinary gauges [35], the rcnormalization of v is not multiplicative in terms of the 
Higgs-boson wave function rcnormalization Zh, but necessitates an additional (divergent) 
shift, Av, i.e., v — Z^{v^ + Av) [36, 37]. Instead, in the HLET-preserving gauges, only the 
Higgs wavefunction Zh is needed to renormalize v, while the shift Av is finite {Av vanishes 
in the MS scheme [38]). Being able to renormahze v only by means of Zh is essential in the 
study of the gauge- (in) dependence of tan/5 in the context of the HLET-preserving gauges, 
appearing in the next section, as well as the definition of effective charges for the Higgs 
sector of the theory [39] . Third, by including fermions in the model, we set up the stage for 
the calculations that will be presented in the next section. 

In section VI we address the issue of the gauge-dependence of the tan /3, one of the most 
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important parameters appearing in multi-Higgs models [40, 41]. We show by means of ex- 
plicit two-loop calculations that the condition of having a vanishing Av does not guarantee 
the GFP-independence of tan /3, because the difference of the corresponding anomalous di- 
mensions displays an explicit dependence on the GFP. This is to be contrasted to what 
happens in the context of the usual (non-HLET-prcscrving) gauges, where the corre- 
sponding difference is GFP- independent, but the GFP-dependence enters again into tan/3 
through the non- vanishing Av [42, 43]. In the calculations presented in this section we make 
extensive use of the £)-formalism, which, together with the HLET, simplifies substantially 
the algebra involved. Our conclusions and outlook are discussed in section VII. Finally, in an 
Appendix we list the Feynman rules for the Abelian Higgs model, together with a collection 
of results needed for several of the calculations appearing in this article. 

II. THE DISPLACEMENT OPERATOR FORMALISM 

In this section we present the derivation of the aforementioned D-formalism. For this 
pourpose it is sufficient to consider a simple scalar theory; the generalization to more 
involved field theories (such as gauge theories) is straightforward, and will be done at the 

end of this section. 

Let then T^n be a bare, one-particle irreducible (IPI) n-point Green function. After 
carrying out the renormalization programme, in o? = 4 — e dimensions, one will have that 

rr^n{X,m^;fi,e) = (/.^rJ„(AR, m^; /x) , (2.1) 

where A, m^, (f) (respectively Ar, m^, 0r) are the bare (respectively renormalized) parameters 
and dynamical field of the theory at hand [44] , 

= 40R, A = ZaAr, m^ = Zm2ml, (2.2) 

and r^„(AR; /i) represents the renormalized n-legs Green function. Note that our definition 
of the bare parameters cj) and A does not include their naive dimensional scalings, namely 
/x"*^/^ and /i.", respectively. Then, we may rewrite the left-hand side (LHS) of (2.1) as 

0%n(A,m2;At,e) = (0r + 50)"r^n(AR + 5A, m| + 5m^; e) , (2.3) 

where the parameter shifts 

50 = (z| - 1)0R, 5X = (Zx - 1)Ar, Sm' = (Z^^ - l)m|, (2.4) 
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quantify tfie difference of the renormalized quantities with respect to the corresponding 
unrenormalized ones in a given renormalization scheme R. Notice that our formulation of 
renormalization treats fields and couplings on equal footing, i.e., as independent fundamental 
parameters of the theory. 

Our next step is to use a Taylor expansion to trade the combinations 0r + dcj), Ar + 5\ 
and m\ + (5m^, for the renormalized parameters 0r, Ar, and mR. To this end, it is natural 
to introduce the differential displacement operator 

in which the shifts are treated as independent parameters, and only at the very end of all the 
manipulations they will be replaced by their actual values in terms of the renormalization 

constants Z^(Ar, m\; e), Za(Ar, m\; e), and Zj„2(Ar, tBr; e). It is then not difficult to derive 
the following (all-order) master equation 

((/)R + 5(t)Y r^n(AR + 5A, m\ + 5m^; e) = ^ 0^ r^"(AR, m^; e) ^ , (2.6) 

or equivalently 

0^rJ„(AR,m^;//) = (e^0^r^n(AR,m^;//,e)), (2.7) 

where (...) means that the CT parameters 50, 5A, and dirP' arc to be set to their actual 
values according to (2.4) after the action of the D operator. The expression (2.7), being 
an all-order result, can now be expanded to any given order in perturbation theory. This 
means that not only V^n should be expanded starting from tree level, but, accordingly, also 
the shifts 50, 5\ and 5m^, together with the displacement operator D (in the latter case the 
expansion starts at the one-loop level, since the tree level shifts are, of course, zero). 

INow, since the shifts are treated as independent parameters, the displacement operator 
at different perturbative orders commute, and one can use the ordinary Taylor expansion 
for the exponentiation of up to three loops, one then gets 

= 1 + + (L»(2) + lz}(i)2j + (^1^(3) + ^{2)£,(i)+ 1^(1)3^ + (2.8) 

with 

L»H = 50(")-^ + 5\^-^^ + 5m2(")^ . (2.9) 
<90R <9Ar dm^ 

The parameter shifts and Sm^^"-^ are loop- wise defined as follows: 

50(") = Z|(")0R, 5A(") = Ar, 5m'^^^ ^ Z^::^ml. (2.10) 



By acting with the operator on the IPI correlation functions 0^r^n(AR, m^; /i, e), we can 
easily determine the expressions for the renormalized correlation functions 0^r^„(AR, m^; /i) 
at the one-, two- and three-loop level, which read 



0^rJW(AR,m^;//) = (i)W0^r(°2(AR,m^;/i)) + 0^rW(AR, m|; e) , 

+ 0Rrin(AR,m|;/x,e). (2.11) 

It is easy to generalize the above formalism to more complicated theories, such as gauge 
theories (with or without SSB) . Let us consider a field theory, where the dynamical degrees 
of freedom, e.g., scalar, spinor or vector fields, are all denoted by (pi, while the static param- 
eters of the theory, such as couplings, masses, VEVs of fields, and GFPs, are collectively 
represented by x. Correspondingly, let us assume that the dynamical and static parameters 
renormalize multiplicatively as follows: 

1 

= ^<|i^iR' ^ = ^xXr. (2.12) 
For such a theory, the displacement operator D will be given by [45] 

^ = E^'^^^ + E'^^^' (2-13) 

^ d(l)iR ^ dxR 
where the CT shifts are defined according to the loop-wise expansion 

^ oo ^ oo 

5<Pi = {Zl-l)<p,n = Sx = {Z,-l)xn = J] ^^'^ ^« • (2-14) 

n=l n=l 

Then, given any IPX Green function Fj-j^^i, it will renormalize according to the formula 



n &n,4.:^(^^-^ 1^) = (e"" n &ii,4>T('^^;i^,e)). (2.15) 

i i 

Equation (2.15) is the master formula of the D-formalism and one of the major results 
of this paper. An immediate consequence of the D-formalism is the renormalization- 
scheme independence of the LHS of (2.15). Specifically, if R and R' denote two arbitrary 



(a) (6) (c) 




id) ie) if) ig) 




(cl) (c2) (c3) (c4) 



FIG. 1: One and two- loop propagator and vertex corrections in the A(/>^ theory. Dots represents 
(j)^ vertices, while crosses and squares represent the one-loop propagator and vertex counterterms 
respectively. Crossing diagrams are not shown. 



renormalization schemes, the renormalization-scheme independence of the bare parameters, 
i.e., (pi = 0jR + S(f)f^^ = 0jR/ + 6(j)f^ ^ and a; = + 6x''^^ = xr' + 6x^^'\ implies that 

n'^"Rrn>?(^R;/^) = H'^rA'r^'^.^r^^R';/^)- (2-i6) 

i i 

Thus, we may employ the Z)-formalism to go from the scheme R to the scheme R', and so 
generally establish the precise relation between the respective IPI Green functions evaluated 
in two different renormalization schemes. 



III. EXAMPLES 

To get acquainted with this new renormalization method, in this section we are going 
to reproduce known results, and clarify the connection between the D operator and the 
conventional renormalization procedure. 
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A. The case of Xcf)'^ 



Let us start by first considering the usual theory defined by 



(3.1) 



At the one-loop level, the first equation of (2.11) gives the usual renormalization prescription 
in the MS scheme, 



— I 



d 

— r^2HAR,m|;/x,e)| 

T^r^4(AR,m|;/x,e) , 



m 



ry(AR,m^;/x,e) 



(3.2) 



R 



where the "bars" indicate that only the infinite part of the corresponding bare Green function 
should be considered. A straightforward calculation of the one-loop diagrams of Fig. 1 (with 
d = 4 — e) gives then 



0, 



z^l 



1 



-A, 



1 



-3A 



(47r)2e ' (47r)2e 
More interesting is the two- loop case, for which the second equation in (2.11) gives 

2 



(3.3) 



2ZI^ 



d 



+ 



i 



D(i)r2 (Ar, ml; e) + (Ar, mi; /x, e) 



D(^)T^'} (Ar, ml; /x, e) + T^'} (Ar, m|; /x, e) 



(1) 7(1) 



m 



R 



D(i)r2(AR, m|; e) + r2)(AR, m|; e) 



(3.4) 



where is now the displacement operator after the one-loop parameters shift has been 
substituted, i.e., 

n(l) _ 7(1) \„ ^ , 7(1)^2 n (1) 



\R ^„.j^ 

with n is the number of external (p legs. Notice that these equations tell us a rather non- 
trivial fact, i.e., that the combination in square brackets on the right-hand side (RHS) is 
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free of overlapping divergences. This is quite striking, since in the calculation of the two 
loop Green function F^^i only IPI bare diagrams must be considered (see Fig. 1 again) and 
no two-loop CT diagram has to be taken into account. Evidently, the D^^^ operator, when 
applied to the lower order (one-loop) F^i\ will generate the necessary CT at the two- loop 
order. 

Since the D operator and the integration over virtual momenta commute, one has two 
possibilities: first carry out the integration and then apply D, or, vice versa, first apply D 
and then integrate. These two approaches have both advantages and disantvantages. In the 
first case, the main advantage is that all explicit references to CTs is removed, and one needs 
to compute IPl irreducible diagrams only (no CT diagrams). The disadvantage is related 
to the fact that, since the D operator involves differentiation with respect to all parameters 
of the theory that undergo renormalization, we need to maintain the explicit dependence 
of the various Green functions on them. This is particularly relevant in the case of gauge 
theories, where one would be faced with the complication of computing diagrams keeping 
the GFP ^ arbitrary [46]. In the second case, instead, one will algebraically generate the 
CTs and recover the conventional formulation. Of course, after applying the D operator, 
one can work at a fixed ^ (say, = 1). To see how the CT diagrams get generated by the 
displacement operator, observe that the one-loop CT to be added to the Lagrangian are 
given by 

'^^i = 77^ ^ ~77^3A2, (3.6) 
(47r)^e (47r)^e 

which are related to the renormalization of the operators 0^ and 0'^, respectively. Then, by 
applying the D operator on the bare one-loop IPI diagrams (a) and (d) of Fig. 1, it is easy 
to see that 

(1) (1) _ A f d'^k 1 1 f d'^k 1 

(cl) + (c2), 

d'^k 1 . f d'^k 1 

{2t:Y (A;2 - m^] 

= (c3) + (c4). (3.7) 



Therefore, irrespectively of the order, one will get the results 



(47r)4l2e' (Att)^ 12e J ' ^ (Att)^ \e'' 12e 



(3.8) 
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B. Massless QED 



We turn to the case of QED, where for simphcity the electron mass is considered to be 
zero to all orders. In this case the one-loop D operator of the photon vacuum polarization, 
assumes the form 

Z,« = Z«e.^.Z-"««|- + Zi". (3.9) 

where e is the QED couphng, and Z^^^ , Z^^^ and Z^^^ are, respectively, the one loop coupling, 
wave-function, and GFP renormalization constants. It is then straightforward to establish 
that the action of the operator D^^^ on the one-loop photon vacuum polarization ^^Ia^ 
vanishes. To begin with, the vacuum polarization is independent of the GFP to all orders. 
In addition, the Abelian gauge symmetry of the theory gives rise to the fundamental Ward 
identity 

q'rl{p,p + q) = S^\p + q)- S^\p) , (3.10) 

where r° and So{k) are the unrenormalized (all-order) photon-electron vertex (with an e 
factored out) and electron propagator, respectively. The requirement that the renormalized 
vertex = -^iF" and the renormalized self-energy S — ZJ^Sq satisfy the same identity im- 
poses the equality Zi — Zf, from which follows immediately that Z^ — Z^^^"^, and therefore, 
after expanding perturbatively, 2Zi^^ + Z^^ — 0. Then, using that, to the given order. 



^^^^'Ia. - 2FW^^ , (3.11) 



we find immediately that D^^^T^^^j^^^ — 0. Evidently, the action of D^^^ produces no two- loop 
CT contributions. This is consistent with the fact that the standard CT graphs (first row 
in Fig. 2) add up to zero, again by virtue of the aforementioned Ward identity (see for 
example [47]). 

At higher orders, the renormalization of the photon vacuum polarization requires CTs 
that originate from two-loop and higher-order self-energy graphs. Likewise, it is not difficult 
to show that, on account of the relation Ze — Z"^^"^ ^ the CTs from one-loop graphs vanish 
identically to all orders. In QED with massive quarks and leptons, mass CTs from all lower- 
order vacuum graphs, including one-loop graphs, contribute to the renormalization of the 
photon vacuum polarization. In this case, the D-formalism is a very practical method to 
reliably calculate such effects. 
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(b) 




FIG. 2: One-loop CTs for massless QED [(a) and (6)], and pure Yang- Mills theory in the back- 
ground field type of gauges [(c) and (d)]. In the latter case we only draw gauge-fixing term insertions 
resulting from the renormalization of the GFP. 



C. Gluon self-energy in the Background Field Method 

We next turn to the case of quark-less QCD, formulated in the BFM, and apply the 
/^-formalism to study the structure of the two-loop CTs appearing in the calculation of the 
background gluon self-energy, to be denoted by T^^^^^. The BFM is a special gauge- fixing 
procedure, which preserves the symmetry of the action under ordinary gauge transformations 
with respect to the background (classical) gauge field A^, while the quantum gauge fields 
appearing in the loops transform homogeneously under the gauge group [48]. As a result, 
the n-point functions (0|r (a;i)A^2 (X2) . . . A^;Jx„) |0) satisfy naive, QED-like Ward 
identities, instead of the Slavnov-Taylor identities valid in the usual covariant quantization. 
Notice also that the n-point functions depend explicitly on the quantum GFP C,q- 

In this case, the one-loop D operator for the (background) gluon vacuum polarization is 
given by 



Here and in the following, stands for the renormalized GFP associated with the quantum 
gauge field A'^{x). Because of the aforementioned background gauge symmetry, Zg and Z-^ 

— 1/2 

satisfy to all orders the QED-like relation Zg = Z^ , from which follows that the first and 
third term of D^^'^ cancel against each other, as happens for the photon vacuum polarization. 
Therefore, D^^^ reduces to 
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On the other hand notice that, contrary to the case of the photon vacuum polarization, the 
background gluon self-energy F^^^^ depends explicitly on ^q. This dependence stems not 
only from the tree-level (quantum) gluon propagator appearing in the loop, 

kaky 



9nv - (1 -Cq)- 



(3.14) 



as happens in the case of the covariant gauges, but in addition from the tree-level vertices 
between a background gluon v4^(p) and two quantum gluons A\^{q) and A\{k), 

^t:M ^) = ar'" [g^p {p-k- ig^q)^ + g.p {k - q)^ + g^,{q-p + ^-^k) ^] , (3.15) 

which, unlike the usual case, also depend on ^q. Finally, as a consequence of the non- 
renormalization of the longitudinal part of the gluon self-energy to all orders, we know that 
= Za [19], where Za is the wave-function renormalization of the quantum gluons; it 
coincides with the one computed in the covariant gauges. Thus, at one-loop (recall that in 
our conventions d = 4 — e) 

(1) q^Ca fl3 \1 

Therefore the action of D^^^ will be non-trivial, and will generate automatically all nec- 
essary two-loop CTs (Fig. 2 second row). Notice that this is in complete accordance with 

the fact that the only CTs appearing in the classic two-loop calculation of [18] are precisely 
related to the gauge-fixing. The advantage of the D-formalism in this respect is that the 
presence of such CTs does not have to be deduced, but is given directly from the very form 
of the D operator. 

In order to verify that the action of D^^^ on F j'* ^ generates indeed the two-loop CT 
contributions, one may follow two equivalent ways. First, one may differentiate directly 
on the Feynman diagram determining F j-* ^ , taking into account the dependence of the 
propagators and vertices on ^g, as shown above, and setting eventually = 1. Thus, one 
generates immediately the Feynman graphs appearing in the second row of Fig. 2. 

The second, more algebraic way, is to actually act with the D^^^ of (3.13) on the expression 
that one obtains for F j'* ^ after the loop integration, and verify that it is actually equal to 
the sum of the CT diagrams, whose individual expressions, calculated at 1, are listed in 
the Table 1 of [18] [where the two diagrams (c) and {d) are called (/) and (m) respectively]. 
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The latter sum reads 

(c) + id) = 1^5"^ (q\. - f ^ . (3.17) 

In doing so, note that, as is well known, the ^Q-dependent part, f ^ > is finite, and is given 
by [49] 

f(l) _ 9^£Axabf^2^ ^ ^ V^Q - 1)(7 + Cq) /0 10^ 

^S^l^^^^ ~ (4^ 1^9 ^/.. - ^/.^^^ j 4 ■ (3-18) 

At this point, acting with D^^^ on f ^ of (3-18), and setting afterwards — 1, we 
immediately recover the result of (3.17), as announced. 



IV. DEFORMATIONS OF THE NIELSEN IDENTITIES 

One of the main advantages of the D-formalism is that it allows to obtain complete 
control over the deformations, caused by renormalization, on the NI, which describe the 
GFP dependence of the bare Green functions. 

To be specific, let us consider the Abehan Higgs-Kibble model in the symmetric phase, 
described by the Lagrangian (see also Appendix A) 

£ = A + £gf + Up . (4.1) 
Here, Ci is the U{1) invariant term 

A = -^F'^^F^^ + (D^$)*(D'^$)-m2$*$-A($*$)^ (4.2) 

with 

F^, = d^A, - d,A^ , D^ = d^-zgYA^, $ = -L(i/ + ^G') , (4.3) 

and the hypercharge assignment — 1. The generic gauge-fixing and Faddeev- Popov ghost 
terms are given by 

£gf = ^B^ + BF, 

£fp = -c{sF) , (4.4) 

where B is an auxiliary non propagating field (which can be eliminated through its equations 
of motion), F is the gauge- fixing function, and s is the BRS operator, giving rise to the 
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following field transformations 



— dfj^c, sH — —gcG, sG — gcH, 

sc = 0, sc = B, sB = 0, (4.5) 

with s nilpotent. Evidently, C is BRS invariant. 

Let us now promote the constant ^ to a (static) field, and introduce an associate anti- 
commuting BRS source 77 (with 77^ = 0) [50] through the transformations 

— T) , ST) — . (4.6) 

The BRS invariance of the original Lagrangian C is then lost, unless we add to it the term 

£n = IvcB , (4.7) 

which couples the 1] source to the other fields. With the Lagrangian term £]<s included, it 
can be easily checked that the so-quantized Lagrangian, 

£ = £1 + £gf + £fp + £n , (4.8) 

is invariant under the extended BRS (eBRS) transformations given by (4.5) and (4.6). If we 
integrate out the auxiliary field B, we find 

B^~F-^r]c, £^^ = -Lf\ c^ = -L^-cF. (4.9) 

In this case, however, the eBRS transformation of the antighost c is not nilpotent, because 
one has 

s^c = I (-sF + ^r]F] , (4.10) 



e V 2e 

where the term in parenthesis is precisely the antighost field equation of motion. Hence, in 
this case, the eBRS algebra closes on-shell. 

The extended Slavnov-Taylor identities one obtains from the eBRS invariant Lagrangian 
are precisely the NI. Specifically, assuming that the B field has been integrated out, the 
generating functional for connected Green functions of our model can be written as 



= J [dA^] [dc] [dc] [dH] [dG] exp |i J d^x ^£ + J'^A^ + J^c + JgC + JhH + JqG 
+Kh(sH) + Kg(sG) + K,(sc)\ } , (4.11) 
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where the eBRS sources K (sometimes called the antifields) have been included only for fields 
transforming non linearly under the algebra of (4.5) and (4.6). Notice that the antighost 
source must be included, since the B field has been integrated out [and therefore c trans- 
forms non- linearly under eBRS, see (4.9)]. Observe that, given a field $, its corresponding 
BRS source obeys opposite statistic, and has ghost charge Qy-p{K^) = — Qfp(^) — 1; 
therefore one has Qfp{Ka^, Kh, Kq, Kc) — {—1, —1, —1, 0} (that Kc has no ghost charge is 
evident from the fact that it couples io sc — B and the action must have zero ghost charge) . 

Thus, eventually the statement of eBRS invariance imphes the following Slavnov-Taylor 
identity: 

Z[J^, Jc, Jc, JiJ, Jg, i^c, Kh, Kq] v] 

= Z[J'^,Jc-ujdf,J'^,Jc,JH,JG,Kc + u;Js,KH + u;JH,Ka + ujJG;^-(jjri,'n], (4-12) 
which, after a Taylor expansion in u; (with o;^ = 0), reduces to the (all-order) NI 

Denoting by <^ all the fields of our model, i.e., Lp = {A^, c, c, H, G}, one can express the 
previous Slavnov-Taylor identities in terms of the effective action ]r[(pci, Kip;^,f]], given by 
the Legendre transform 

T[ip^i,K^;i,r]] = Z[J^, K^;i,r]] - J d'^xJ^f^i , (4.14) 

with (fci — 5Z/SJ^. Taking into account the relations 

^r-'- '^--'^'^ 

we obtain 

ST ^ ST ST ST ST ST ST ^ 
J^y^^ + + + + ^^^^ = ° • ^'-''^ 

Finally, after differentiating with respect to 77 and setting 77 afterwards to zero (with 77^ = 0), 
we obtain the NI in its final form 



Introducing the so-called Slavnov-Taylor operator Sr, the NI (4.17) is often written down 
as 



= Sr 

r]=Q 
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dr,T 



(4.18) 



-9«r^= + 

t ^ 




FIG. 3: Graphic representation of the NI (in this case the one for the tadpole). Notice that the 
arrows indicate the flow of the ghost charge. The last term is absent to all orders in the gauge. 

where we have used that [S^ , 9^] = 0. It is now easy to see that, due to ghost charge conser- 
vation, the first term on the RHS of (4.17) becomes relevant only when one is calculating the 
gauge-dependence of ghost Green functions. For instance, it contributes a term — S^Fg^^^ 
for the NI of the ghost self-energy Fee. 

A graphical representation of the NI can be obtained as follows. The t] coupling can 
be represented by a heavy dot attached to a dotted line with an arrow that indicates the 
flow of the ghost charge. Likewise, the coupling of the BRS sources Kh and Kq should 
by represented by a (S> symbol attached to an arrowed dotted line. On the other hand, the 
coupling of Kc should be represented by ® attached to a dashed line without an arrow, since 

has no ghost charge. In this way, a consistent graphical picture conveying information of 
the flow of the ghost charge emerges. For example, the NI for the Higgs tadpole, that reads 
(assuming CP conservation) 

—d^^H = —St 

can be represented as in Fig. 3. It can be explicitly checked that the last term on the RHS 
of (4.19) is zero to all orders in perturbation theory in the conventional gauge- fixing 
scheme. 

As already mentioned in the Introduction, unlike the Slavnov-Taylor identities which 
are derived from the BRS invariance of the theory, the NI are obtained from the eBRS 
symmetry; therefore they will not remain unmodified by the process of renormalization, but 
they will be deformed. In what follows, we are going to determine a closed expression for 
the deformation of the NI under renormalization, by applying the Z)-formalism. 

Let us therefore denote by x the parameters of the theory (coupling, masses, etc.) except 



rjc-L KsH 



(4.19) 



18 



of the GFP ^, and with 0j the fields A^,!!, and G, with renormahzation relations given 
by (2.12). To elucidate our points, we will initially consider Green functions involving n 
fields of one type, say 0, and then we will generalize our results to Green functions involving 
different types of fields. Based on the basic formulae (2.6) and (2.7) of the D-formalism, we 
can now study the response of the renormalized Green functions to a variation A^r of the 
renormalized GFP ^r, i.e., 



0RrJ.(a;R,eR + AeR;/x)-0Sr?.(a;R,eR;/x) = cj>lA^Kd^^T^n{xK,^R; l^) + 0(Aei) • (4.20) 



In stating (4.20), we have assumed that the renormalized fields and parameters 0r and 
are GFP-independent in a given renormahzation scheme R. In fact, this is the case, if the 
renormalized parameters are evaluated from renormahzation conditions that are manifestly 
GFP independent, e.g., from physical observables, such as 5'- matrix elements, or from gauge- 
invariant operators, within a good regularization scheme that preserves the Slavnov-Taylor 
identities. In general, the GFP independence of the renormalized parameters in a specific 
scheme can be determined in two ways: (i) by showing that the corresponding CTs are GFP- 
independent; (ii) by comparing it to another gauge-invariant scheme, such as the MS [38] or 
the pole-mass scheme [51, 52]. An exception to the above are the renormalized field 0r and 
its VEV vr. These parameters are GFP-independent, but their respective renormalizations 



and are in general GFP- dependent, even within the MS scheme. Thus, although we 
will initially assume that 0r and xr are GFP-independent, we will also comment on the 
modifications that should be considered, if xr were depending on ^r. 

Let us concentrate on the LHS of (4.20), and more precisely on its first term. From (2.1), 
we will have that 



Substituting for the shifts their expressions in terms of renormahzation constants [cf. (2.14)], 




+ 




(4.21) 
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and then inserting the result back into (4.21), we obtain for the RHS of (4.20) that 



<pR9^R^^n{xR,CK;ii) = 0" a^j,(Z^^R) d^r^n{x,^;ii,e) + Z^cp'^ld^^Z^^ XB^d^r^n{x,C; ii,e) 

+ 710^ (9felnz|) rJ„(xR,eR;//) , (4.22) 

which is equivalent to the identity 

— V "I — / \ 

^&rJn(2;R,CR;/") = Z^^d^j^{Z^^R)^d^T^n{x,^;ii,€) + [d^^Z^jxRd^T^n{x,^; n,e) 

+ n (d^. In 4) rJ„(a:R, ^r; pt) . (4.23) 

Finally, including ^ in the parameters x, we may cast (4.23) into the slightly more compact 
form: 

n — — / \ 

9^K^<l,n{x^;n) = Z^Z^d^r^n{x;ii,e) + Z^ i^d^^Z^jxRd^r^n{x;ii,e) 



+ n{d^^lnZ-?)TUxR;f^). 



(4.24) 



A similar line of arguments can be followed, if xr depends on ^r. In this case, the factor 
{d^^Z,j:)xji contained in the second term on the RHS of (4.24) should be replaced with 
55j,(Z^xr), for xr ^ ^R. 

We may now rely on the D-formalism to express the RHS of (4.24) entirely in terms of 
renormalized parameters. In this way, we obtain the deformed NI 

n_ / r 

9^R^^n (xr; h) = Z 2 Z^{e Sr 9^r<^» (xr; r], e 



(4.25) 

+ Z^ (^d^^Z^^ xr ^e'°9j;j^r^n(xR;/x, e) ^ + n [d^^lnzj) rJ„(xR;/.). 

Notice that, up to an overall constant, only the first term on the RHS of (4.25) can be related 
to the undeformed, unrenormahzed NI (4.18), where the IPI Green functions involved are 
evaluated with renormalized parameters. The appearance of the other terms is a consequence 
of the renormalization process. We emphasize that, unlike earlier formal treatments [26], 
Equation (4.25) furnishes for the first time the precise closed and, in principle, calculable 
expression for the deformations of the NI. 

It is interesting to remark that even the action of the operator present in the first term 
on the RHS of (4.25) gives rise to further deformations. To make this last point more explicit, 
we first observe that the action of the D operator in (4.25) is actually equivalent to the action 
of a reduced operator D, without the field derivatives, i.e., D — D — S4>i d/dcpiR. Making 
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use of this last observation, we decompose e^, or equivalently e^, as = 1 + (e^ — 1). 
Subsequently, noticing that the operators D and St commute, i.e., \D ,8^] — 0, the first 

term on the RHS of (4.25) can be rewritten, up to the multiplicative factors Z^, as 



e ot 



9^r^n(xR;/i,?7,e) 



(4.26) 



Thus, only the first term can be expressed in terms of the undeformed, unrenormalized Nl 
of (4.18), whereas the second one is an additive deformation of the NI that results in from 
a BRS variation of another function [26]. 

The generahzation of (4.25) to Green functions involving different types of fields is 
straightforward [53], and reads 



^^^^Yii " n ^'Pi ^« \ [ dr, Tj^. (xr; /X, 77, e) 

i 

i X 

+ J2ni(d^JnZl^T^^^^,{xn;^) ■ 

i 

In order to get a feci on the structure of Equation (4.27), let us apply it to the lowest non- 
trivial order. Expanding consistently, according to the D-formalism, we obtain the following 
one-loop result: 

i 

i X 

+ E (^&^^'^) x^d,,T^^]^.^{x^) , (4.28) 

X 

where the first term is the one- loop NI of (4.17), with the simple replacement of bare 
parameters by renormalizcd ones. 

A further simplification of the formulae above occurs in a gauge-invariant renormaliza- 
tion scheme, such as the MS scheme. In this case, all terms proportional to d^j^Z^, for 
which X is related to a gauge-invariant operator {e.g., a gauge-coupling constant or a gauge- 
invariant mass parameter), will drop out from the RHS of (4.27), as their multiplicative 
renormalization constants will be ^R-independent. 
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V. THE ABELIAN HIGGS MODEL IN THE HLET-PRESERVING GAUGES 



In this section we will discuss the Abelian Higgs model quantized in the type of gauges 
that preserve the Higgs-boson Low Energy Theorem (HLET). The Lagrangian defining the 
model is given by 

£ = A + >Cgf + ^Cfp . (5.1) 
Here Ci is the gauge-invariant part of the Lagrangian 

A = -^F^^F^^ + (Z}^$) -m2$*$- A ($*$)^ 

i i 

- V2hJ,^^f^ - x/2/i2/2^$7i' - \/2/ii/«$7^ - V2h2fJ^f^ . (5.2) 

In (5.2), F^j, = d^Ai, — d^^A^ is the U{1)y Field strength tensor, = d^ — igYA^ is the 
corresponding covariant derivative, and 

$ = -^[v + H + iG^ (5.3) 

is the complex Abehan Higgs field, composed from the two real fields H and G, where 
($) = v/\/2 is its VEV that signifies the spontaneous symmetry breaking of U{1)y- The 

hypercharge quantum numbers of the different fields are assigned according to Y$ = 1, 
y/' = —Y2 = 1, and = Y^ = 0. For simplicity, we finally assume that the Yukawa 
couplings hi are real. 

For the gauge-fixing term Cqf in (5.1), we will adopt the one introduced in [27, 34], i.e., 

Cgf = -^[d,A'^ + ^glm{<^>')r 



a 



d,A^ + ^g{v + H)G \ (5.4) 



which in turn induces the Faddeev-Popov ghost term 

£pp = -c [a^ + 2/^Re($2)] c 

= -c{92 + /e[(^ + ^)'-G'^]}c. (5.5) 

Notice the presence of the two GFPs, ^ and a. The gauge-fixing scheme described by £gf 
and the BRS- induced ghost term L^v (which will be referred to as the scheme), in ad- 
dition to the technical advantages mentioned in [27], belong to the special class of schemes 
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which respect the so-called HLET [28-32] for off-shell unrenormahzed Green functions be- 
yond the tree level [34] (see also our discussion below). The reason is that the complete 
Lagrangian, including the gauge-fixing and ghost sectors, is invariant under the transla- 
tional transformation [34]: 

H ^ H + a, V ^ V - a. (5.6) 
For example, the conventional gauge, described by the gauge-fixing term. 



r - 1 



1 



2 



d.A'^ + -^^gvlm{^) , (5.7) 



violates this translational symmetry (5.6) due to its explicit dependence on v. As a conse- 
quence of this violation, the HLET given below by (5.10) is no longer valid beyond the tree 
level. 

The translational symmetry (5.6) is sufficient to show the vahdity of the HLET to all 
orders. As a result of this symmetry, the entire effective action T[H, v] (where we suppress 
all other fields and parameters) satisfies the identity 

T[H,v] = T[H + a,v -a], (5.8) 

which immediately implies the translational Ward identity 

5T dT 

In general, upon n functional differentiations with respect to the field we obtain the 
HLET for off-shell (unrenormahzed) n-point IPI Green functions Vh^: 

Th^.. = ^ . (5.10) 

In the above formula, the Higgs-boson insertion in Fj^n+i is evaluated at zero momentum. 
This result should be contrasted with the one obtained in the i?^ gauge, in which the 
HLET described by the relation (5.10) will be grossly violated by gauge-mediated quantum 
effects [34]. 

The translational identity (5.9) has an immediate consequence on the way the Higgs VEV 
is renormalized. The most general way of renormalizing v is given by [36, 37] 

V = Z^v^ = Zl{v^ + ^v) . (5.11) 
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It should be remembered that Av differs from 5v defined previously in (2.14), since Sv — 
(Zy—l) v^. The quantity Av may be split into two parts: one part that contains the divergent 
contribution proportional to 1/e, to be denoted by Af |div, and a finite, renormalization- 

scheme dependent piece, to be denoted by At>|fin, i.e., Av = Af |div + Av\fi^. The crucial 
point is that if the HLET is exact, then one must have that Af l^iv = 0. To prove this, we 
observe that the Higgs tadpole and the effective potential F, by virtue of (5.9), satisfy 
the equality Fh — dyV. Given that = Z^Fh, we find 

F| = Z^d^T = ziz-'d,^T = z|z-^a,,F^ , (5.12) 

where in the last step wc have used the fact that F(a;; /x, e) = F^(a;R; /x). Therefore we get 
the condition 

Ziz-^ = Cfin , (5.13) 

where Cfin is an UV finite constant. In perturbation theory, this finite constant may be 

decomposed as Cfin = 1 + ^cgn, where 5c^n is a higher-order scheme dependence. In the MS 

scheme, the multiplicative renormalization constants have no UV finite pieces, so that one 

1 

would get unavoidably cgn = 1, with 5cfai — 0, and therefore Z^ — Z^, which is equivalent 

1 

to Av|div = to all orders, with v — Z^vyi- In other renormalization schemes, one needs to 

impose that (5.10) holds true after renormalization, so that again cgn = 1; this can be done 

1 

without intrinsic inconsistencies in the HLET preserving gauges, by imposing that Zfj — Zy, 
even for the UV finite pieces. 

In the remainder of the section, we will present the Lagrangian of the Abelian Higgs 
model, after the SSB of U{1)y- This will enable us to set up the stage for the two-loop 
calculations related to the issue of gauge dependence of tan which will be discussed in the 
next section. The full Lagrangian of the U{1)y model after SSB may be written down as a 
sum of four terms: 

f ■'■ \ 4 ■'-2 2 

Lq — — Xv m V , 

4 2 

£i = -{Xv^ + m^v)H , 

^ = lid.Hf - ^{3Xv' + m')H' + ^-{d.Gf - ^[(A + a^g'y + m']G' 
- \f,,F^^ - ^^{d.A^f + \g'v''A,A^ -{a- l)gvG{d,A^) 
+ a^ca^c - ig^v^cc + fii^Yd,, - hv)U , 

i 
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+ g'vHA^A^^ + ^g'H'A^A^ + \g^G^A,A^ + (a + l)gGA,{d^H) 

+ (a - l)gHA^{d^G) - 2Cg%Hcc - Cg'H'cc + ^g'G'cc 

+ gA^ Yl fi^h^^PLfi - J2 ^^f^Hf, - J2 yl^hJil.GU , (5.14) 

i i i 

where all the fields and parameters are bare. Observe that a = 1 is the only value that avoids 
the appearance of mixed propagators between the would-be Goldstone boson G and the 
gauge boson A^. We will therefore renormalize the model by imposing the latter condition 
on the renormalized GFP itr, i.e., — 1. This condition becomes rather subtle at two loops, 
and especially when the D operator is applied after integration, in which case one needs to 
keep the (one loop) full dependence on a in all the quantities under study. However, as we 
already mentioned, our calculational task may considerably be simplified if the D-formalism 
is applied before the integration of the loop momenta. 

In the HLET-preserving gauges, the propagators take on the following form: 







i 






ml 


i 





G 



k"^ — m\ 



-Qt^y +1 



^ \ k^ki^ 



a J k"^ — -ml 

(7 C 



^i^k) = , (5.15) 

p-mi 

where the particle masses are related to the independent parameters of the theory through 

= g^v^ , = 3Av^ + m? , = (A + (7^g^)v^ + m^ , 

"^c — ) — hiV . (5.16) 

The complete set of Feynman rules is finally given in Fig. 6 of the Appendix. 



VI. GAUGE DEPENDENCE OF TAN /3 AT TWO LOOPS 

In models with two elementary Higgs bosons, Hi and H2, one of the fundamental pa- 
rameters is the ratio of their vacuum expectation values, vi and U2, respectively [40, 41]. In 
particular, the quantity usually denoted as tan/?, is defined at tree-level as tan/? = v^jvx. 
When quantum corrections are included tan/3 develops a non-trivial dependence on the 
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renormalization mass /i, as well as on the unphysical GFP. Given that tan P is exten- 
sively used in parametrizing new physics effects in many popular extensions of the Stan- 
dard Model, such as two-Higgs models and almost all supersymmetric versions, this type of 
gauge-dependence is an undesirable feature. Various studies have therefore addressed the 
question under which conditions a gauge-independent definition of the running tan (3 could 
become possible [42, 43]. 

In studying these issues, there appears to be a subtle interplay between being able to 
set A.v\iiiy — and showing that the difference ^yni — TiJa of the anomalous dimensions 
is independent of the GFP. In this section we will explore in detail this connection, and 
demonstrate that, contrary to what one might naively have expected, it is not possible to 
establish the gauge- independence of — 1h2-i least not within a conventional field- 
theoretic framework. 

The basic observation which suggests a link between the gauge-independence of tan/3 
and Ai>|div = is the following. If Ai>i|div = 0, for i = 1, 2, then vi — Z^.Viji, and therefore 



tan /3 renormalizes as 



tan/3 ^ ^ = ^ ^ = ^ tan/3R . (6.1) 



The renormalization group equation for tan /3r is thus given by 

d tan /3r 



dt 



= (iH^ - IH^ tan/^R, (6.2) 



where t = In^u, and 7//^ = — | J^(ln Z^:/.) is the anomalous dimension of the Higgs field Hi. 
If at this point one could show that, in the class of gauges where Afj|div = 0, the difference 
(7i?2 ~ 7-ffi) is GFP-independent, one would have a solution to the problem. The crucial 
point in this argument is precisely that the two conditions need be satisfied simultaneously. 
Indeed, having a gauge-fixing scheme where {'yH2 ~1hi) is GFP-independent does not imply, 
by virtue of (6.2), the GFP-independence of tan/^R, unless one could demonstrate that, 
within the same scheme, one is also able to set Afj|div = 0; if the latter condition cannot 
be enforced the renormalization group equation that tan/^R satisfies is simply not that of 
(6.2), since the starting assumption is not valid. Satisfying both conditions simultaneously 
is far from trivial. In fact, as we will see in detail, in the context of both the and 
the HLET-preserving gauges, these two conditions cannot be simultaneously met, for 
different reasons. Specifically, in the R^ gauges each of the two-loop anomalous dimension 
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7ifj, consists of two pieces: (i) a gauge dependent 0{g'^) polynomial, common to both, and 
(ii) a gauge independent 0{h^g^), which is different for Hi and H2. Thus, in taking the 
difference {jh2 ~1hi), one finds a GFP-independent answer for this quantity [42]. However, 
since the gauges are not of the HLET-preserving type, one cannot set Afj|div = 0, and 
therefore (6.2) receives additional (gauge- dependent) contributions. On the other hand, the 

gauge preserves the HLET by construction, and one may set Avj|div — 0, thus enforcing 
the validity of (6.2); however, as we will see in detail in what follows, the two-loop calculation 
reveals that, in these latter gauges, (7^/2 — TifJ is in fact GFP-dependent. 

In order to demonstrate this, it is actually sufficient to consider the Abehan-Higgs model, 
despite the fact that it contains only one Higgs field. The rationale is that in the context of 
-R^ gauges the 0{h'^g'^) contributing to the Higgs boson anomalous dimension turns out to 
be GFP-dependent. Therefore, given that, in general, each Higgs boson couples differently to 
the fermions {i.e., the corresponding Yukawa couplings are independent parameters), even 
if there were a second Higgs boson, this gauge dependence could not in general cancel in the 
difference {'^h^ ~ IHi)- 

In the rest of this section we will prove the gauge dependence of the 0{h?g^) contributions 
to 7)^ . To that end, we will employ two independent, but complementary, approaches. In 
the first approach we will exploit the validity of the HLET in the R^ gauges in order to 
eventually obtain the (non- vanishing) first derivative of Z^^ with respect to ^ from the two- 
loop effective potential, without actually computing Higgs-boson self-energies. Second, we 

will explicitly compute Z\j' from the corresponding two-loop diagrams. Since for determining 

(2) 

Z)j one needs to consider only the divergent terms proportional to the external momentum 
g^, whereas terms whose dimensionality is saturated by masses do not contribute to Z^^\ 
one may carry out the calculation in the symmetric phase, when v — Q. In both approaches 
we will employ the D-formalism developed in the previous sections in order to enforce the 
necessary cancellations of the overlapping divergences, without explicit reference to CT 
diagrams. 

Let us first write down how the relevant Green functions renormalize. Prom the second 
equation of (2.11), one finds the following relations 

i {ZHZm^m'i + Z]^Zx\rvIY'^Kb. = |^^^^riJ^(a;R; e) r2^(a;R; e) } 

I. J Ictr,=1 

i{ZHZ^.ml + 3ZlZxXnvlY^^ = \Dmr^^l{xn;pi,€) + r^^l{xK;pi,e)]\ 
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6i 

6i 



l(TR = l 



{ZjjZxfhnVn = {L'Wr^J) (xr; e) + rg^l^R; /x, e)} | 

(Z|fZ;,)^'^AR = |D(i)r«(xR;/x,e) + rgl(xR;//,e)|| , (6.3) 

I. J l(TR=l 



where xr collectively denotes all the renormalized parameters of the model, and, according 
to our definitions, the one-loop displacement operator is given by 

r)(i) _ ^7(1),, ^ I 7(i)\ ^ ^ 7«^2 9 yd) 9 r,ii). d yd) d 



d n 



+ (6.4) 



with n the number of external Higgs legs. We emphasize that the above equations are 
not independent from each other, since, due to the HLET, they are related by successive 
differentiation with respect to v. Notice also that, on the RHS of (6.3), all overlapping diver- 
gences should cancel, a fact which furnishes a very stringent check of the entire calculation. 
From the above equations one can infer the gauge-dependence of the Higgs wave-function 
renormalization constant at two loops. Specifically, expanding the last equation of (6.3), we 
find 



2Z1J) + Zf = -(Z«)^-2Z«Z«-^{Da)rg(xR;/.,6) + rg(xR;/.,6)}[ . (6.5) 



Differentiation with respect to ^, taking into account that, in the MS scheme that we are 
using, Zx is GFP independent, yields the following identity 

= (4" + + ~ {fl<'>r«(x„;M,.) + rg(xH;^.)}[_^^_ . (6.6) 

One may easily verify that similar equations can also be obtained starting from any of the 
first three equations of (6.3). 

Due to the HLET, from the diagrams contributing to the effective potential we can 
extract information about the Higgs tadpole, mass, tri- and quadri-linear couplings, by 
simply differentiating with respect to v. Moreover, since we are interested only in the 
contributions which depend on the Yukawa couplings, we only need to consider the two-loop 
fermionic effective potential contributions, shown in Fig. 4. Introducing the integrals 

rf^fc f (ft 1 



{2'kY J {2'kY (A;2 - ml){P - m|)[(A; + ^)2 - ml] 



n=i,j,k 



J' 

+ 
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(a) (6) 



FIG. 4: The two loop fermionic contribution to the effective potential. Due to the HLET, differ- 
entiation with respect to the Higges VEV v provides the one-, two-, three- and four-Higgs Green 
functions (at zero momentum). 



.2e 



(2^ J (27r)^(P-mf)(£2-m2; 



l + ^(2-ln 



+ . . 



(6.7) 



where In/z = ln(47r/i ) — 7^;, and the dots stand for finite parts, one finds 



--hjd 
2 ' 



^Hf, - ^Lfj^ + (2ml - ]^m]j 1 K 



-i^d{d-2) 



(6.8) 



The contributions to Tj^i which depend on the Yukawa couphngs are simply obtained by 
differentiating the above expressions four times with respect to v. As far as the term D(^)r^4 
is concerned, one can use the results of Appendix A for the one-loop effective potential and 
renormalization constants (notice that in this case, one cannot limit one's attention to 
the fermionic contributions only, due to the dependence on the Yukawa couplings of the 
renormalization constants). 

The combination of these two terms leads to a massive cancellation, yielding finally 

29' 



(47r)4e 



(6.9) 



(2') 

thus establishing the GFP-dependence of Z\j'^^ in the gauges. 

As a check for the consistency of the procedure, we can evaluate the full two-loop Higgs 
wave-function renormalization constant of our model, through the direct calculation 
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FIG. 5: Fermion diagrams contributing /i? terms to the two-loop Higgs wave-function renor- 

(2) 

malization constant in the symmetric phase. 

of the relevant Feynman diagrams of the Higgs-boson self-energy, shown in Fig. 5. As 
mentioned earlier, we will calculate in the symmetric phase, f = 0, keeping only contributions 
proportional to the Yukawa couplings hi 

We will use directly the D-formalism to explicitly check that all overlapping divergences 
cancel, and to get the fermionic contributions to Z\j' through the formula 

^^'^ = {Wr«(a:K;/.,6) + Fg(a:K;/.,6)}[ . (6.10) 



Once again one should keep in mind that F^a must be calculated at a general value of a 
(which means that the (a — 1) part of the A^GH vertex will also contribute). The final 
result is given by 

which coincides with (6.9) after differentiating with respect to ^. From this expression one 
can determine the two-loop anomalous dimension of the Higgs field, which is given by [54] 

d 



(2) 

n 



where x„ denotes collectively the free parameters of the theory, p„e their mass-dimension, 
and C]j the coefficient of the simple pole of the Higgs IPI self-energy. Therefore 



(2) 



^ d Id 1 V-.. ^ ^ ^1 



Evidently, despite the fact that Av = exactly, the two-loop running of tan (3 turns out to 
be GFP-dependent. 
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VII. CONCLUSIONS 



We have developed a new formalism for determining the renormalization and the GFP- 
dependence of Green functions to all orders in perturbation theory. The formalism makes 
use of the fact that the renormahzed Green functions are obtained by displacing both the 
unrenormalized fields and fundamental parameters of the theory with respect to the renor- 
mahzed ones. Due to this property, wc have called it the displacement operator formalism 
or, in short, the D-formalism. With the help of this formalism the CTs necessary for the 
renormalization of any Green function can be unambiguously determined, to any given or- 
der of perturbation theory. In particular, if one applies the D-operator before integrating 
over the loop momenta, one can systematically generate all the CTs that would have been 
obtained in the conventional diagrammatic framework. We explicitly demonstrate the full 
potential of the /^-formalism by considering several known examples of renormalization of 
theories up to 2- loops, such as a 0^-theory, QED, and QCD in the BFM gauge. 

One of the great advantages of the D-formalism is that it can be used to calculate the 
precise form of deformation of symmetries, which are modified in the process of renormal- 
ization, such as the NI. Hence, the dependence of the renormahzed Green functions on the 
renormalized GFPs can be computed exactly, thereby offering a new method for evaluating 
the GFP-dependcnce within a given scheme of renormalization, e.g.^ the MS scheme, the OS 
scheme, etc. Given that a concrete, closed formula describing the deformation of the NI to 
all orders exists now, one should be able to conclusively settle various formal issues related 
to gauge-invariance. Most notably, it would be interesting to revisit the important question 
of the all-order gauge-invariance of the pole of the unstable particles, together with other 
topics related to the gauge- invariant formulation of resonant transition amplitudes [55] . 

In theories with SSB, in addition to the Slavnov-Taylor identities and the NI, the HLET 
plays an important role as well. The ordinary gauge violates the HLET for off-shell 
IPI correlation functions. In order to explore the constraints imposed by the HLET we 
have resorted to a toy field theory with SSB, the Abelian Higgs model, which was quantized 
using a HLET-preserving gauge. An important consequence of these gauges is that the VEV 
V of a Higgs field renormalizes multiplicatively by the Higgs wavef unction, so there is no 
additional UV infinite shift to f , i.e. Afdiv = 0. Employing the D-formalism in the context 
of a HLET-preserving gauge, we have shown that the fundamental quantity tan /3, defined 
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in the two-Higgs models, is GFP-dependent at two loops, exactly as happens for the usual 
gauge, in which there exist additional UV infinite shifts to the Higgs VEVs. The analysis 
presented here strongly suggests that the tan /3 cannot be made gauge-independent within 
the framework of conventional Green functions, even if the gauge-fixing employed respects 
all relevant symmetries and constraints, such as the HLET. These results motivate one 
to explore the possibility of defining tan/? at higher orders through the GFP- independent 
effective Green functions constructed by means of the pinch technique [56]. In particular, 
it would be interesting to extend the concept and construction of the Higgs-boson effective 
charge [39] to the case of multi-Higgs models, and more specifically to supersymmetric 
theories. We hope to be able to report progress on this subjects in the near future. 

The formulation developed in this article presents novel perspectives for the study of 
several other known topics. Specifically, the D-formalism may be used to systematically 
investigate the renormalization-scheme dependence of IPI correlation functions. It may 
also be employed to algebraically determine the restoring terms of a "bad" UV regularizing 
scheme, i.e., a scheme that does not preserve the Slavnov-Taylor identities. Since it provides 
all-order information on the renormalization of Green functions under study, it might be 
useful in controlling the calculation of non-perturbative efi'ects, such as those related to the 
dynamics of renormalons [57]. 

The D-formalism can be straightforwardly extended to systematize the procedure of 
renormalizing non-renormalizable field theories. In particular, it may be used to organize 
the infinite series of CTs needed to renormalize such theories. But even in the case of 
renormalizable perturbative field theories the D-formalism can be automated, for example 
with the aid of a computational package, to reliably compute all the CTs required for the 
renormalization of IPI correlation functions at high orders. It would therefore be interesting 
to explore these new horizons opening up, embarking into a study of the onset of non- 
perturbative dynamics at very high orders of perturbation theory [58] . 
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APPENDIX A: FEYNMAN RULES AND RENORMALIZATION OF THE 
ABELIAN HIGGS MODEL 



The Feynman rules derived from the Lagrangian density given in (5.14) are hsted in 
Fig. 6. The model can be renormalized with the renormalization condition (Tr = 1, by 
introducing the following renormalization constants 



= Zm?m\, A = ZaAr, g = Zgg^, ^ = 



V = Z^Vr, 



(Al) 



Notice that left and right fermions get renormalized with different renormalization constants. 

The determination of the renormalization constant above, is simplified, due to the HLET: 
for example, the knowledge of the fermion self-energy will automatically imply the knowledge 
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FIG. 6: The Feynman rules for the Abelian Higgs model in the HLET type of gauge. 
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i { (lOA^ + 2£,g'^\ + - 4^,. hf) + (8A + 2£,g'^) m'^v'^ + 2m^] 



^ { (lOA^ + 2e.9^A + 3ff4 - 4 hj) + (4A + ^g^) m^v} 



Th 



^{[(6 + 205^-4E,/^?] 

-(<SA + 2C'r)///- - (OOA- + 12Ai;//- + 18y^ - 24 5], ''"} 
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+ (8A - 6^.g2)TO2 + (20A2 + AX^g^ + 6g* + 6(^g* -8J2^ K) ^'^l 
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TABLE I: The one-loop divergent parts of the Green functions of the Abelian Higgs model in the 
HLET preserving gauge (with a = 1). 

of the Higgs-fermion-fermion vertex through the differentiation of the former with respect 
to the Higgs vev v (notice that, in this particular case, we would find immediately that this 
vertex is one- loop finite within the HLET gauges) . 

In table I, we report all the divergent parts for the one-loop Green functions of the model. 
Using these results, one finds 

(6 + 2ey-4E.^^ 



h'Zh = 1 + 
fi'Zo = 1 + 
H'Za = 1 + 
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(47r)2e 




1 + ^ 

(47r)2e ' 
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(47r)2e 
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(47r)2e 




4A+(_6 + 6ey + 4E>? 
(47r)2e 



A l^i "'I 



(A2) 



Notice that Z\Zg = 1 as it should due to the Ward identities of the theory. Combining the 
divergent parts of the Green functions given in table I, together with the above renormal- 
ization constants, one can explicitly check the validity of the renormalized NI of Eq.(4.28) 
at one loop. 

The effective potential is particularly useful in the context of the HLET-preserving gauges, 
furnishing a substantial amount of information with relatively moderate effort [27]. In 
dimensional regularization, the one-loop effective potential reads 

where In/i^ = ln(47r/i^) — ^e- As stated earlier it is important to keep explicitly the depen- 
dence on (T, since this will play a crucial role in cancelling the overlapping divergences in 
our two-loop expressions, when the D-formalism is applied after integration. 

At two loops, there are two basic topologies contributing to the effective potential, shown 
in Fig. 7. It turns out that all these two loop diagrams can be expressed in terms of the 
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FIG. 7: Topologies contributing to the two loop effective potential, i, j, and k labels all the possible 
field combinations allowed by the Feynman rules. 



integrals Lij and Kijk introduced in Eq.(6.7). In table 2 we report the results of the (scalar 
and vector) diagrams. Using these results in conjunction with Eq.(6.6), one can obtain the 
full gauge-dependence of the two-loop Higgs wave-function renormalization and anomalous 
dimension. Specifically, 

d_ i2) ^ -29' E. - 4A^^ - + 2g' 8Xg' + 24^/ 
^ (47r)4e (47r)4e2 ' 

9 . (2) _ 2g' + 4A^^ + Q^g' - 2g' 
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+ "'Vi"'[iAA + ice ~ 2Lac] + 27^(1 - e)(3™i - 2ml + ^l) [KhAc - Khcc] 

+ ''"\rnf"^' i^HAA + Khcc - 2KhAc] } 
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-iS.g'^Lnc 


{v2)gc 
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TABLE II: Scalar and vector two loop contributions to the effective potential (with a = 1). The 
fermionic contributions are given in Eq.(6.8) of the text. 
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We end by reporting for completeness the full two-loop Higgs wave function renormalization 
and anomalous dimension: 
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